********************************************************************************
* Turnout
********************************************************************************

use dataset_voters, clear

gen post = 1 if year == 2012
replace post = 0 if year == 2000
gen anyfrackingXpost = anyfracking*post
tabulate year, gen(year_)

gen highfracking2 = highfracking
replace highfracking2 = 1 if inlist(di_state, "CO", "WY")
gen highfracking3 = highfracking2
replace highfracking3 = 1 if inlist(di_state, "NM", "OH")

eststo clear
preserve
	keep if highfracking2 == 1
	eststo: reg turnout lnwells year_2 year_3 year_4 year_5, ///
		a(zipcode) cl(zipcode)
	estadd local ZFE = "Yes"
	estadd local YFE = "Yes"
	eststo: reg turnout anywells year_2 year_3 year_4 year_5, ///
		a(zipcode) cl(zipcode)
	estadd local ZFE = "Yes"
	estadd local YFE = "Yes"
	eststo: reg turnout anyfrackingXpost year_4, a(zipcode) cl(zipcode)
	estadd local ZFE = "Yes"
	estadd local YFE = "Yes"
restore

keep if highfracking3 == 1
eststo: reg turnout lnwells year_2 year_3 year_4 year_5, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"
estadd local YFE = "Yes"
eststo: reg turnout anywells year_2 year_3 year_4 year_5, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"
estadd local YFE = "Yes"
eststo: reg turnout anyfrackingXpost year_4, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"
estadd local YFE = "Yes"

label variable lnwells "Log wells"
label variable anywells "Any wells"
label variable anyfrackingXpost "Ever fracked X post"

#delimit;

esttab 
	using "_output/tableA6a.tex"
	,
		replace
		booktabs
		order(lnwells anywells anyfrackingXpost)
		drop(*year* _cons)
		cells(b(fmt(3) star) se(par))
		compress
		label
		gap
		collabels(none)
		nomtitles
		width(\linewidth)
		mgroups(
			"+ WY, CO" "+ WY, CO, OH, NM", 
			pattern(1 0 0 1 0 0)
			prefix(\multicolumn{@span}{c}{) suffix(})
			span erepeat(\cmidrule(lr){@span})
		)	
		stats(
			YFE
			ZFE
			N_clust
			N
			,
				label(
					"Year fixed effects"
					"Zip code fixed effects" 
					"Clusters"
					"Observations"
				)
				fmt(0 %20.0fc)
		)
		;
	
#delimit cr

********************************************************************************
* Contributions
********************************************************************************

use dataset_donors, clear

gen post = 1 if year == 2012
replace post = 0 if year == 2000
gen anyfrackingXpost = anyfracking*post

tabulate year, gen(year_)

gen highfracking2 = highfracking
replace highfracking2 = 1 if inlist(STATE, "CO", "WY")
gen highfracking3 = highfracking2
replace highfracking3 = 1 if inlist(STATE, "NM", "OH")

eststo clear
preserve
	keep if highfracking2 == 1
	eststo: areg ltotamount lnwells year_2 year_3 year_4 year_5 ///
		year_6 year_7 year_8, a(zipcode) cl(zipcode)
	estadd local ZFE = "Yes"
	estadd local YFE = "Yes"
	eststo: areg ltotamount anywells year_2 year_3 year_4 year_5 ///
		year_6 year_7 year_8, a(zipcode) cl(zipcode)
	estadd local ZFE = "Yes"
	estadd local YFE = "Yes"
	eststo: areg ltotamount anyfrackingXpost year_7, a(zipcode) cl(zipcode)
	estadd local ZFE = "Yes"
	estadd local YFE = "Yes"
restore

keep if highfracking3 == 1
eststo: areg ltotamount lnwells year_2 year_3 year_4 year_5 ///
	year_6 year_7 year_8, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"
estadd local YFE = "Yes"
eststo: areg ltotamount anywells year_2 year_3 year_4 year_5 ///
	year_6 year_7 year_8, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"
estadd local YFE = "Yes"
eststo: areg ltotamount anyfrackingXpost year_7, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"
estadd local YFE = "Yes"

label variable lnwells "Log wells"
label variable anywells "Any wells"
label variable anyfrackingXpost "Ever fracked X post"

#delimit;

esttab 
	using "_output/tableA6b.tex"
	,
		replace
		booktabs
		order(lnwells anywells anyfrackingXpost)
		cells(b(fmt(3) star) se(par))
		drop(*year* _cons*)
		compress
		label
		gap
		width(\linewidth)
		collabels(none)
		nomtitles
		mgroups(
		    "+ WY, CO" "+ WY, CO, OH, NM", 
		    pattern(1 0 0 1 0 0)
		    prefix(\multicolumn{@span}{c}{) suffix(})
		    span erepeat(\cmidrule(lr){@span})
		)	
		stats(
			YFE
			ZFE
			N_clust
			N,
				label(
					"Year fixed effects"
					"Zip code fixed effects"
					"Clusters"
					"Observations"
				)
				fmt(0 %20.0fc)
		)	
		star(* 0.05 ** 0.01 *** 0.001)
		;
		
#delimit cr

********************************************************************************
* End
********************************************************************************
